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Abstract 

The effects of surface polar plionons on electronic transport properties of monolayer graphene 
are studied by using a Monte Carlo simulation. Specifically, the low-field electron mobility and 
saturation velocity are examined for different substrates (SiC, Si02, and Hf02) in comparison 
to the intrinsic case. While the results show that the low-field mobility can be substantially 
reduced by the introduction of surface polar phonon scattering, corresponding degradation of the 
saturation velocity is not observed for all three substrates at room temperature. It is also found 
that surface polar phonons can influence graphene electrical resistivity even at low temperature, 
leading potentially to inaccurate estimation of the acoustic phonon deformation potential constant. 

PACS numbers: 72.80.Vp,72.10.Di,72.20.Ht,73.50.Dn 
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Since the experimental realization of graphene in 2004,- its extraordinary electric proper- 
ties have led to near exponentially growing attnetion in possible applications. In particular, 
the linear band structure in the vicinity of the Dirac points, and resultant massless fermions, 
have excited interest in both its applicability for electronic devices and as a testbed for 
quantum field theory phenomena.- The extremely high intrinsic mobilities that have been 
observed in suspended graphene at room temperature,- in conjunction with its true two di- 
mensional structure of graphene, naturally lead to its employment in higher speed graphene 
transistors devoid of short channel effects.^ For such high frequency applications, the sat- 
uration current and low-field mobility are material properties that need to be accurately 
determined. 

When a graphene layer is in close proximity to a polar substrate, inelastic carrier scat- 
tering with surface polar phonons (SPPs) can result in significant reduction in the low-field 
mobility of graphene." Due to the inelastic nature of SPP, they also provide pathway to cur- 
rent saturation, in conjunction with, or as a substitute for, intrinsic optical phonons. While 
a number of studies generally suggested the negative influence of reduced saturation veloc- 
ities (due to the relative small energies of SPPs),^"^^ conflicting reports exist that predict 
a very different picture (including enhanced velocities) based on the analysis of Boltzmann 
transport equation.-ii^ 

In this work, we investigate the influence of carrier scattering due to SPPs on the elec- 
tronic transport properties of monolayer graphene, via a full-band ensemble Monte Carlo 
method. Particularly, the low-field electron mobility and saturation velocity are calculated 
in the presence of three different substrates (SiC, Si02, and Hf02) and compared with 
those of intrinsic graphene at room temperature. Furthermore, we examine the impact of 
SPPs on the low temperature electrical resistivity, with attention to its implication on the 
experimental determination of the acoustic phonon deformation potential constant. 

The Monte Carlo model adopted in the calculation utilizes the complete electron and 
phonon spectra in the first Brillouin zone. While a tight-binding band is used for the elec- 
tronic energy structure, all six branches of the graphene phonon spectra are considered with 
the phonon dispersion and electron-phonon scattering rates obtained from first principles 
calculations.— The effect of degeneracy is accounted for by the rejection technique, after final 
state selection.— The distribution function is obtained self- consistently from the ensemble 
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simulation. The SPP scattering rate is introduced by following^i^ 



Tsiki) h ^ 2e{qy q 

where q = |kf — ki| is the SPP momentum, Ef (Ei) the final (initial) electron energies, d 
the distance between monolayer graphene and the substrate (0.4 nm), us the SPP energy, 
and J-"^ = — -i^jr^) the square of Frohlich coupling constant. The (1 + cos 9) factor 

originates from the overlap integral of the pseudospin part of the electron wavefunction and 
e{q) is the dielectric function in graphene. In addition, J-" contains the dependence on high 
(low) frequency dielectric constant of the polar substrate (k^.) along with the normalized 
area A. The scattering by ionized impurities in the substrate is not included in the effort to 
clearly identify the role of SPPs. The values for relevant substrate parameters can be found 



in Refs. 



6- 



Figure 1 shows the SPP scattering rates calculated at 300 K for three substrates, SiC, 
Si02, and Hf02. For comparison, the strength of electron-optical phonon interaction inher- 
ent in graphene (via the deformation potential) is also plotted from a first-principles analy- 
sis.— The results clearly illustrate the dominance of SPPs over intrinsic optical phonons in 
graphene for the entire electron energy under consideration. Furthermore, this enhancement 
in scattering is more pronounced for the substrates with small us that is apparent from the 
onset energy of emission process (e.g., Hf02 - 19.4 meV; Si02 - 60.0 meV; SiC - 116 meV). 
Due to the Coulombic nature, the SPP scattering is a function of electron density n and sub- 
sequent screening in graphene. Throughout the calculation, we assume n = 1 x 10^^ cm^^ 
along with the static screening function e{q) in the random-phase approximation. 15 The 
intrinsic scattering via deformation potential interaction has no dependence on n. 

Figure 2 shows the electron drift velocity vs. the electric field at T=300 K, for intrinsic 
graphene and graphene on different substrates. As expected, the addition of electron-SPP 
interactions lead to general decrease of electron drift velocities in the low-field region. Conse- 
quently, the low-field mobilities are reduced for all three substrate with the largest decrease 
in Hf02. However, the difference with the intrinsic case becomes progressively smaller with 
the increasing field (thus, the increasing average electron energy) and, in the case of Si02 
and SiC, the velocity appears to saturate at a higher value. No degradation is observed even 
for HfOo. Clearly, this behavior does not follow the saturation velocity model suggested 
by Ref. |9| based on the onset of optical phonon emission, intrinsic phonon or otherwise 



(i.e., SPPs). Rather, it is in general agreement with more recent studies that solved the 
Boltzmann transport equation numerically assuming a displaced Fermi-Dirac distribution 
function.-ii^ 

The apparent inconsistency of increased total scattering rate and absence of velocity 
degradation in the high-field region may be explained by considering the linear energy dis- 
persion and the characteristics of the Frohlich interactions. As it is well known, a change 
in electron energy (say, via an inelastic scattering) does not automatically relax the drift 
velocity in monolayer graphene near the K and K' points. What matters is the direction 
of the final momentum. Due to the Coulombic nature, the Frohlich interactions, including 
those by SPPs, prefer small angle events. When a SPP scattering occurs, the electron will 
likely emit (absorb) a SPP, with the phonon momentum pointing along the radial direction 
toward (away from) the Dirac point. Consequently, the velocity of the scattered electron 
may not change significantly. In contrast, optical phonon scattering of intrinsic graphene is 
via deformation potential interactions and, thus, randomize/relax effectively the direction 
of the electron momentum/drift velocity. Ignoring the difference between the Frohlich and 
deformation potential interactions^ can lead to inaccurate depiction of transport properties. 

When the electron-SPP interactions are taken into account along with other scattering 
mechanisms, they provide an additional channel for energy and momentum relaxation. Ac- 
cordingly, the average electron energy becomes substantially lower at a given electric field. 
As the applied field increases, it means that the shift in the distribution function (along 
the direction of the electric field) would be smaller in energy along with a shorter tail com- 
pared to that in intrinsic graphene. This is clearly visible in Fig. 3 plotted for the ky = 
cross-section at 20 kV/cm. The observed shorter tails in the negative space (i.e., with 
negative velocities) as well as in the high energy region with nonlinear dispersion^ appear 
to more than compensate the additional momentum relaxation of SPP scattering on the 
Si02 and SiC substrates, leading to larger saturation velocities shown in Fig. 2. In the case 
of Hf02, the distribution is much less heated with a significant population in the negative 
half space. Due to the very strong inelastic scattering by surface phonons, the drift velocity 
may not have reached the saturation point at 20 kV/cm and could grow further. Cooling 
of electron distributions via surface phonon interactions may also reduce the self-heating as 
the generated heat (i.e., SPPs) is on the substrate and, thus, can be more readily dissipated. 
Another interesting point to note from Fig. 3 is that the electron distribution function in 
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graphene does not resemble that of a highly degenerate case, certainly not when subject 
to an appreciable electric field. A simple approximation of a displaced box-like function^^ii 
cannot describe transport properties accurately. 

The strong infiuence of the substrate on low-field mobility observed at 300 K raises a 
possibility that the SPP scattering may be an efficient mechanism even at low temperatures. 
This is a distinct possibility due to the small phonon energies us- Figure 4 plots the electrical 
resistivity p in graphene vs. temperature for all four cases under consideration. Indeed, the 
resistivity is substantially affected by the substrate conditions at T > 50 K and the effect 
of SPP scattering can not be immediately separated from those of other mechanisms (most 
importantly, the acoustic phonon scattering). The acoustic deformation potential is often 
determined experimentally utilizing the assumption that the low temperature resistivity 
is linearly dependent on temperature with a slope that is proportional to the square of 
acoustic deformation potentiaL-ii^ However, the calculation results indicate that this may 
not be possible; the extracted deformation potential (Dae) is not an intrinsic quantity but an 
"effective" parameter that includes substrate effects and screening (i.e., the charge density 
in graphene) among others. 

Concerning the specific results, the intrinsic case gives Dac ~ 6.8 eV with a well-defined 
linear region. This (i.e., the strong linearity) is due to the negligible contribution of optical 
phonon scattering at low temperatures with relatively large energies {huop ~ 160 meV). 
The small deviation of Dac from the deformation potential constant obtained directly from 
the acoustic phonon scattering rate (~ 4.5 eV)-^ may be attributed to a degenerate electron 
density considered in the present calculation (n = 1 x 10^^ cm^^). It is also interesting to note 
that this value (6.8 eV) is actually rather close to the experimentally extracted on a non-polar 
substrate (7.8 eV).— For the SiC substrate, we estimate Dac ~ 7.1 eV similar to the intrinsic 
result. As T ^ 150 K, however, there is a substantial increase in SPP scattering and the 
difference between the SiC and intrinsic cases becomes more discernible. When Si02 is used 
as the substrate, the slope of p is further increased and becomes nonlinear earlier due to the 
small SPP energy. The deduced value in the linear region (T ^ 125 K) gives the appearance 
of Dac ~ 13.2 eV, which is not unlike 16—18 eV estimated experimentally on Si02."^ Finally, 
IIf02 has the largest effect among the three substrates. The resistivity is greatly increased 
with no apparent linear region in the temperature range under consideration. Accordingly, 
it is not plausible to determine Dac- The relevance of SPP scattering even in the T ^ 150 K 
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regime can explain, at least in part, the very disparate results for the magnitude of acoustic 
phonon deformation potential in graphene.-'^'i^'i^ 

In summary, we investigate the effects of substrate-induced SPP scattering on the elec- 
tronic transport properties in monolayer graphene by using a full-band ensemble Monte 
Carlo simulation. It is found that the electron velocity-field characteristics are highly de- 
pendent on the choice of substrate at room temperature. Specifically, the Frohlich nature 
of the interaction appears to be crucial for correctly describing the saturation of drift ve- 
locity. The simulation also shows that the SPP scattering remains an efficient mechanism 
even at low temperatures (T ^ 50 K), substantially affecting the electrical resistivity. This 
makes it difficult to experimentally determine the acoustic deformation potential for intrinsic 
graphene in close proximity to a substrate. 

This work was supported, in part, by the Defense Advanced Research Projects 
Agency/HRL Laboratories (the CERA program), the US Army Research Offce, and the 
SRC Focus Center on Functional Engineered Nano Architectonics (FENA). MBN wishes 
to acknowledge partial support from the Offce of Basic Energy Sciences, US DOE at Oak 
Ridge National Lab under contract DE-AC05-00OR22725 with UT-Battelle, LLC. JMZ ac- 
knowledges support from NSF under the IR/D program. 



6 



^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grig- 

orieva, and A. A. Firsov, Science 306, 666 (2004). 
^ See, for example, A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. 

Geim, Rev. Mod. Phys. 81, 109 (2009). 
^ K. I. Bolotin, K. J. Sikes, Z. Jiang, M. Klima, G. Fudenberg, J. Hone, P. Kim, and H. L. 

Stormer, Solid State Commun. 146, 351 (2008). 

F. Schwierz, Nat. Nanotechnol. 5, 487 (2010). 

5 J. Chen, C. Jang, S. Xiao, M. Ishigami, and M. S. Fuhrer, Nat. Nanotechnol. 3, 206 (2008). 

6 S. Fratini and F. Guinea, Phys. Rev. B 77, 195415 (2008). 

^ A. Konar, T. Fang, and D. Jena, Phys. Rev. B 82, 115452 (2010). 
8 V. Perebeinos and Ph. Avouris, Phys. Rev. B 81, 195442 (2010). 

^ I. Meric, M. Y. Han, A. F. Young, B. Ozyilmaz, P. Kim, and K. L. Shepard, Nat. Nanotechnol. 
3, 654 (2008). 

° J. Chauhan and J. Guo, Appl. Phys. Lett. 95, 023120 (2009). 

^ A. Barreiro, M. Lazzeri, J. Moser, F. Mauri, and A. Bachtold, Phys. Rev. Lett. 103, 076601 

(2009). 

2 A. M. Dasilva, K. Zou, J. K. Jain, and J. Zhu, Phys. Rev. Lett. 104, 236601 (2010). 

3 K. M. Borysenko, J. T. Mullen, E. A. Barry, S. Paul, Y. G. Semenov, J. M. Zavada, M. 
Buongiorno Nardelh, and K. W. Kim, Phys. Rev. B 81, 121412(R) (2010). 

^ X. Li, E. A. Barry, J. M. Zavada, M. Buongiorno Nardelh, and K. W. Kim, Appl. Phys. Lett. 

97, 082101 (2010). 
^ E. H. Hwang and S. Das Sarma, Phys. Rev. B 75, 205418 (2007). 
6 E. H. Hwang and S. Das Sarma, Phys. Rev. B 77, 115449 (2008). 

^ X. Hong, A. Posadas, K. Zou, C. H. Ahn, and J. Zhu, Phys. Rev. Lett. 102, 136808 (2009). 



7 



Figure Captions 



Figure 1. (Color online) Surface polar phonon scattering rate of graphene electron on SiC, 
Si02, and Hf02, for the electron density n of 1 x 10^^ cm~^ at 300 K. Also plotted is 
the intrinsic graphene optical phonon scattering rate at 300 K ("intrinsic") obtained from 
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Figure 2. (Color online) Electron drift velocity in graphene on SiC, Si02, and Hf02. The 
intrinsic case without substrate is also shown. The electron density is 1 x 10^^ cm^^ at 
300 K. 

Figure 3. (Color online) Electron distribution function with different substrate conditions: 
intrinsic (dash-dotted), SiC (solid), Si02 (dashed), and Hf02 (dotted), for the ky = cross- 
section at 20 kV/cm. The box-like function corresponds to the Fermi-Dirac distribution 
displaced by the Si02 SPP energy (60 meV) in a simple, metal- like approximation (i.e., Ep 3> 
ksT) with n = 1 X 10^^ cm^^. For comparison, the equilibrium Fermi-Dirac distribution at 
300 K is also plotted. 

Figure 4. (Color online) Electrical resistivity vs. temperature with different substrate 
conditions: intrinsic (circle), SiC (square), Si02 (triangle), and Hf02 (diamond), with 
n = 1 X 10^^ cm~^. The slope of the straight lines are used to extract the acoustic de- 
formation potential. 
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FIG. 3: Li et al. 
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